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SUMMARY 

The operation of a shore-based automatic digitally 
recording water quality monitor on the St, Clair River down- 
stream from Sarnia during the period April 1971 - July 1972 
is discussed. This monitor measures chloride, temperature, 
dissolved oxygen, pH and conductivity. Instruments of this 
nature employed singly are limited to detecting extreme 
values of chemical parameters which may be missed during 
routine sampling, as well as interparameter correlations and 
periodic variations. As data are collected at a fixed point 
some other checking is required to determine if the point is 
representative of the river quality. The data translation, 
conditioning and analysis requirements are extremely demanding 
on both personnel and computer time for such limited information 
Operational problems with the intake and shore mounted equip- 
ment, particularly the mechanical components caused numerous 
breakdowns or resulted in unusuable data. 

Typical ranges of data obtained v/ere as follows: 
chloride, 20-4 mg/1 ,• dissolved oxygen, 7 0-110% saturation; 
pH, 7,5-9.0 units; and conductivity, 200-350 ymhos/cm. 

High values of chloride data collected during April 
1971 correlated with meteorological data, and strongly 
suggested runoff of salt from city streets. In April-June 
1:971, the results obtained were characterized by generally 
non-normal frequency distributions and diurnal variations of 



temperature, dissolved oxygen and pH consistent with solar 
heating and photosynthesis. In November 1971-January 1972, 
higher probabilities of normal distribution accompanied by 
the almost complete absence of periodic variations indicated 
the interaction of a large number O'f phenomena, including 
possibly both natural processes and industrial loadings. 
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WATER CHIMISTIY DATA FOR THE ST. CLAIR 
RIVER AT CORUNITA, ONTARIO',, AS 
DETERMINED BY CONTINOOUS WATER QUALITY MONITOR 

1. INTRODUCTION 

Determining water quality by grab sampling tech- 
niques has inherent problems associated with defining a 
significance of the sample. Normally, numerous. s.amples are 
collected to evolve a statistical representation. However, a 
valid statistical value for bodies of water which are gener|ally 
both spatially and temporally variant may not be obtained from 
a few random samples. Nor are estimates o^f extremities likely 
to be available from such data. 

Automatic water quality monitoring equipment which 
could measure and record various water quality parameters at 
fiKed intervals was immediately recognized as being valuable 
in defining statistical and temporal relationships. Equip- 
ment of this type was used on the Potomac and Ohio rivers in 
the early 19i0's (McDermO'tt, 1969; Klein, Dunsmore and Horton, 
1968) for' water management purposes. However, to date, tech- 
nology has been unable tO' provide the instr'umentation required 
tO' measure many important water quality parameters. Nevertheless, 
automatic instr'uments can be used tO' record some of the simpler 
to measure parameters. It has thus been possible to^ record 
means,, standard deviations, dietributio'ns or periodicities with 
respect to time, as well as extremities which would likely be 
missed during routine sampling. Even this in,forma,tion has been 
difficult to obtain , as is evidenced by the volume O'f equipment 
sales compared tO' the minute quantities o^f published data avail- 
able. There is no doubt that the data are mandatory at least on 
a limited scale to^ design surveys and develop' both statistical 

and mathematical water managero>ent models, 
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Submersible water quality recording meters, have 
been operated in lakes Erie and Ontario (Palmer and Izatt, 
1970; Palmer, 1972). These units have the advantages of 
portability and easy calibration and electronic repair in 
the laboratory. On the other hand, miniaturization results 
in a loss of reliability in the sutaiersible meter. Calibra- 
tions cannot be adjusted in situ. Also, malfunction of any 
probe is not readily detected during the period of operation. 
A shore-based meter overcomes these difficulties, and in 
addition allows cleaning of the probes during operation and 
rapid calibration checking without significant down-time. 

The St. Clair River was selected as a location for 
trial of a shore-based meter as it is influenced by the heavy 
density of chemical industries in the Sarnia area- To this 
end, a continuous water quality monitor with digital recording 
system (robot monitor) was operated from April 1971 to July 
197 2 at Corunna, downstream from the industrial area (Figure 
II, 



2. DESCRIPTION OF EQUIPMENT 

The automatic water quality measuring system used 
in this study was manufactured by Schneider Instruments Co., 
who used the term "robot monitor" to describe the system 
(Model No. RM 25), The monitor measures chloride, temperature, 
dissolved oxyen (DO), pH and conductivity, and records these 
parameters periodically at a set sampling frequency on computer 
compatible punched paper tape . Although the sampling interval 
was generally 1 hour, the system could be set for several 
intervals between 10 minutes and 2 hours. 

Temperature is measured by a silicon resistance 
element, which is incorporated in a DC bridge. The conductivity 
of the water between two platinum electrodes is measured with 
a 60 cycle AC bridge. Automatic temperature compensation 
produces a result corrected to 25°C. pH is measured with a 
glass electrode and silver-silver chloride reference electrode 
using a KCl solution as a salt bridge, and chloride with a 
similar system containing a chloride-sensitive glass electrode . 
DO is measured with a polarographic electrode, which consists 
of a gold cathode and four silver anodes immersed in a KCl 
solution, separated from the solution being analyzed by a 
membrane permeable to oxygen. Molecular oxygen passes from 
the water through the membrane, and is reduced at the cathode 
by electrons from the external voltage supply. These electrons 
are returned to the circuit by silver dissolution at the anodes. 



resulting in a current which is directly proportional to 
the number of oxygen molecules passing through the membrane. 
During the period of this study, the robot monitor 
was housed at the Corunna sewage treatment plant. The intake 
was initially located at a depth of 1.5 m, approximately 10 m 
from shore, but as this was blocked by eKtensive weed growth 
in July 1971 and not very representative of the river, it 
was relocated about 30 m from shore at a depth of 3.5 m. 
A detailed description of the installation, alo^ng with problems 
encountered during the operatio«,, is given, in Appendix I. 



3. STATISTICAL ANALYSIS OF DATA 

Stochastic methods like time eeries and Markov 
chain analysis are valuable in proceeeing large amounts of 

data produced by the robot monitor. Time series analysis 
detects periodic variations related to solar heating and 
photosynthesis, lake , seiches and periodic waste loadings. 
A frequency and time distribution analysis defines the 
frequency of occurrence of extreme data values which may 
be detrimental to water use or aquatic life. Markov chain 
analysis can show the persistence O'f extreme values O'f para- 
meters, i.e. will an extreme value remain present for a long 
period of time? 

Prior to the application of these methods,, data 
preparation is necessary. As the records, contain errors 
derived from improper tape punching or false signals from 
the instrumentatio'n o^r senso'rs, data editing is performed. 
For meaningful time series analysis, long term trends must 
be removed by trend correcting, in order that the mean value 
for each parameter may be applicable throughout the whole 
time period, 

3.1 Tape Translation 

The punched paper tape produced by the digitizer 
is read and the data transcribed onto magnetic tape. Dis- 
solved oxygen figures are converted to percent saturation,, 
using saturation data given by Standard Methods (APHA, 196 5) . 
Conductivity figures are standardized to 18:*C using the Great 
Lakes composite conductivity-temperature data of Rodgers (1962) 
The data thus obtained are termed "raw data". 



3.2 Data Editing 

FrO'm the raw data, the mean and standard deviation 
of each chemical parameter (i.e. conductivity, pH, etc) are 
deterinined . All observatiO'ns gxeater o^r less than twO' Btan= 
dard deviations from the mean are rejected and replaced by 
linear interpolatiO'n between adjacent valid readings. This 
criterion allowed the rejection of nearly all the spurious 
data* 

3.3 Trend Correcting 

The entire edited record is divided into a number 
of equal sub-sets. These were approximate integer multiples 
of a day, in order to preserve diurnal variations (Palmer 
and Izatt, 1970'). A polynomial is fitted to the means of 
these eubsets against time , with the maximum order being two 
less than the number of subsets, or 5, whichever is less, A 
lower order is used if no improvement of fit occurs. A random 
scatter of the division means is poorly represented by poly- 
nomial fitting; consequently, the polynomials obtained are 
subjected tO' an F-test for significance of fit. If the pody- 
nomial is significant at 5 percent,, a predicted value for the 
chemical parameter is obtained for each data point and sub- 
tracted from the reading. The difference is added to the 
grand mean of all the intervals to obtain a value which illus- 
strates variations whose periode are less than the period of 
the subset. 



3.4 Frequency and Time Distribution Analysis 

The range of values obtained in each set of raw, 
edited and trend corrected data is divided into ten equal 
segments, and the frequency of observations falling within 
each segment determined. The frequency area moments are calcu- 
lated. The distribution indicates the probability of 
obtaining high or low results, while the moments provide a 
measure of the mean, variance, skewness and kurtosis of the 
distribution (Snedecor and Cochran, 1967). 

As many natural processes tend to produce a distri- 
bution which is normally or log-normally distributed with time, 
•a measure of the probability of finding such distributions is 
desired. Four methods of calculating such a probability were 
tested. These include (1) the Kolmogorov-Smirnov test and 
(2) the Von Mises test applied to the continuous distribution 
of data points; as well as (3) the Von Mises discrete test 
and (4) the Chi-squared test applied to the frequency table 
already obtained. Data obtained during this comparison are 
presented in Appendix 3. As the result obtained for a theoreti- 
cal normal distribution by the Von Mises continuous test differed 
from that expected by the authors, and the discrete tests gave 
variable results depending upon the manner in which frequencies 
were expressed, the Kolmogorov-Smirnov test was selected for 
use in time distribution analysis. 

Both frequency and time distribution analysis are 
useful in assessing the magnitude of change obtained during 
data editing and trend correction. 



3.5 Time Series Analysis 

A time series analysis was conducted on the raw, 
edited and trend corrected data for each parameter to obtain 
the autocorrelation coefficients and the smoothed spectral 
density estimates. Analysis at each processing stage allows 
a comparison of raw, edited and trend corrected data to be made, 
and thus provides an estimate of the effects of each step 
on the data . 

The time variation of the autocorrelation coefficient 
indicates the frequency of changes occurring in water chemistry 
parameters, and aids in the interpretation of such changes. 
Spectral analysis improves the resolution of complex periodi- 
cities and is useful when the auto -correlation function shows 
several interfering frequencies. 

For meaningful time series analysis, data must be 
"pre-whitened" or smoothed to remove irregular short-term fluctua- 
tions. Each smoothed data point is obtained by forming a running 
mean of several consecutive observations. Although an equally 
weighted running mean is simplest to calculate, in order to 
obtain a smoothly varying frequency response, the observations 
are weighted proportionally to the binomial coefficients, 
resulting in a binomially weighted running mean (Panofeky and 
Brier, 1968, p. 147-151). In general, hourly data were pre- 
whitened over 5 time intervals allowing the observation of 
spectral peaks at all periods of interest. 
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Autocorrelation coefficients r^^ of the time series 

JLj 

formed by each of the parameters are computed for lag time L 
up to 100 hours with the equation: 



N _ _ 

r^ = Z (X. - X) (X. .^ - X) , . . (1) 

L . n 1 X + ij 

1=1 __________ 

S ^N 



where X. is the ith reading of any parameter 
X is the mean O'f N readings O'f X. 

S ^ is the variance of N readings of X. 

X 1 

Spectral density estimates are obtained by Fourier 
cosine transformation of the autocorrelation function (Panofsky 
and Brier, 1968, p. 143). The spectral estimates B. are 
obtained from.: 



m-1 

B. = k fr + r cos (iiT) + 2E r^. cos(Li7T/m)l ... (2) 
1 — '^ o m - T L / ' J 

m L=l 



where r , r^ and r are autocorrelation coefficients 
o L m 

for lags 0, L and m (m - maximum lag time). 

L and i are lag times, each varying individually 
from to m. 

k is a constant, k = 1/2 for i = and i = m 

k = 1 for other values of i. 



As equation (2) produces rapidly fluctuating spectral 
density estimates, smoothed spectral estimates are computed by 
a technique called Hanning (Blackman, 1958), according to: 

S. ^ 0.25 B. , + 0.50 B. + 0.25 B.^, ... (3) 
1 1-1 1 1+1 
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where S. = smoO'thed spectral estimate. 
These estimates are plotted against frequency ^ j" 

where At is the measuring interval 

(in most cases, 1 hr) and i and m are as in (2). 

In order to determine the statistical significance 
of peaks in the spectral estimate curve, an average line was 
drawn through the points, and the ratios O'f maxima and minima 
to this average were calculated. The limiting levels O'f 

significance for these maxima and minima may then be esti- 
mated from the plot given by Jenkins (1968, p. 82), where the 
number of degrees of freedom v is estimated from, v = (2N - m/2)/m 
where N and m are the total number of readings and maximum 
lag time, as stated previously.. 



3.6 Markov Chain Analysis 

The range between the minimum and maximum readings 
for each chemical parameter is divided into ten states as in 
Section 3.4, and these are numbered consecutively from one to 
ten, O'ne being the state containing the lowest readings. A 
10' X 10 one hO'Ur transition probability matrix T.__ (1) (where 
m = initial state and n = state after 1 hour) is determined 
for the period of record. This matrix defines the probability 
of occurrence of each of the ten pO'Ssible states after one hour, 
given the present state. From this matrix, it is possible to 
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determine the final state pxobability vector after any 
number (k.) of hours, ass'uming that the initial state 
probability vector for k hours is equal to the final state 

vector for (k-1) hours. By comparing the probability after 
k hours with the final state probability vector, one may 
estimate the time required for the system to return to an 
average state following the recording of an extreme value. 

The final state probability vector is the proba- 
bility of each state occurring over the period of the record, 
as determined by frequency analysis (Sec. 3.4) . Probability 
vectors for k hours may be determined by the transition 
matrix using the relation: 

m=i 
where 

q (k) = final state probability vector 



'n 



after k hours 



q (k-1) = probability vector after k-1 
™ hours or initial state prob- 

ability vector for k = 1 

T,_.. 11) = one hour probability transitio^n 
"^ matrix 

The initial state is usually defined at an extreme 
reading (e.g. state 1, probability = 1.0), and probability 
vectors for 1-4 0' hours are determined for each chemical para- 
meter. A rapid progression to the final state frequency vector 

inidcates that the extreme values are rare and unlikely to 
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persist for long periO"ds of tiBT , while a slow progression 
indicates the values vaxy greatly and that extremes are 
not unlikely. 
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4. RESULTS AND DISCUSSION 

The robot monitor was operated continuously from 

April 1971 to July 1972. Readings were made at hourly 

intervals^ except in April 1971 when a two hour interval 
was used. 

In many cases, ^ equipment malfunctioning resulted 
in invalid data. These included the action of weeds and ice, 
as well as calibration difficulties and digitizer failures. 
In August-September 1971, the paper tape was punched off 
center, and could no't be translated. In May-July 1972, 
increasing electronic difficulties in the digitizer resulted 
in bad data being obtained. These problems are discuseed in 
detail in Appendix 1. A complete record of all data tapes 
obtained is given in Appendix 2. 

The only periods in which reasO'nable results were 
obtained were April-June 1971 and November 1971-February 1972, 
except for a period in late November, when several sensors 
were out of calibration. Other time periods will not be 
discussed further. 

4.1 Data Preparation 

Complete listings of raw, edited and trend corrected 
data are rO'Utinely printed so that the effect of each proces- 
sing step on the data record may be ascertained. Table 1 
illustrates a sample listing of data from April 1971, 
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Tables 2-4 and Figures 2-4 illustrate the frequency 
and time distribution analyses for raw, edited and trend 

corrected chloride data from April 1971. The removal of 
miepunched data on the paper tape by editing is illustrated. 
As such mispunching produces either very high or very low 
results, editing often produces a conBiderable change. In 
this case, the range of chloride values is seen to decrease 
from 20-325 mg/1 for raw data, to 20-58 mg/1 for edited data. 
This decrease is also accompanied by a decrease in the momenta 
of the distribution,, as shown in Tables 2 and 3. 

The results of polynomial fitting to the edited data 

fO'r purpo'ses O'f trend correction are given in Table 5. Typical 
plots of polynomials obtained during trend correGtion have been 
published previously (Palmer and Izatt, 1970). In several 
cases, the number of subsets was varied in order to determine 
whether or not any bias occurred in the data from use of a 
given number of subsets. This factor is illustrated in Table 
5 for December 1971. The results are similar for either 10 
subsets (approximately 3 days each) or 14 subsets (approximately 
2 days each) , although increasing the number O'f subsets 
tends to increase the level of significance. 

The effect of trend correcting on the frequency 
and time distribution analysis of chloride data from April 
1971 is illustrated in Tables 3-4 and Figures 3--4. Trend 
correcting results in a further decrease in the range of 
data to 22-42 mg/1. While the practice of trend correction 
is necessary for time series analysis, it should be approached 
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with cautiO'n. This decrease in range, together with removal 
O'f skewnesB and increased probability of normal distribution,, 
actually results in the submergence of extreme values of 
data, which are important ecologically. An example of this 

is the chloride data discussed in Section 4.4. 

Means and standard deviations of all data thus 
obtained are presented in Table 6. Standard deviations 
obtained by robot monitor for chloride, DO and conductivity 
are similar to or slightly lower than standard deviations of 

six grab Bamples obtained between May and November ? 1971. 
As a, fractiO'n of the mean, standard deviations of the warm 

and cold weather data are similar, except for conductivity, 
where they are higher during the cold weather period. 

4.2 Frequency and Time Distribution Analysis 

Results of the Kolmogorov-Smiriiov test for normal, 
log-normal or exponential distribution are presented in Table 

7. Distribution is an impO'rtant characteristic since the 
application of parametric statistical tests to the data assumes 
normal distributiO'n. If a large number of processes, natural 
or artificial, are combining to produce effects ^ a normal or 
log-normal distribution is expected. For exaiiple, conductivity 
is influenced by all ionic species present in the water, which 
in turn may be affected by a considerable number of factors 
such as municipal and industrial discharges, rainfall, biological 
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activity, wave and current actiO'ii. Non=normally distributed 
properties, o^n the other hand,^ frequently result when one or 
two factors are primarily responsible for observed variations 

In general |. cold weather data appear tO' be more 
normally distributed than warm weather data,, particularly in 
the case of co^nductivity and dissolved O'xygen, which are 
stro^ngly affected by photosynthesis in the growing season. 
Distributions for' chloride and temperature were unpredictable 
and the pH distribution was generally no^t normal. 



4.3 Time Series Analysis 

The results O'f time series analysis, presented in 
Table 8, reflect the difference between the warm and codd 
weather records, with warm weather data indicating diurnal 
periodicities attributable to solar heating and photosynthesis 
and cold weather data generally being less periodic in natures 
Typical autocorrelatlO'n functions from June 1971 are presented 
in Figures 5-8 and the variance density spectra derived from 
these in Figures 9-12. 

4.3.1 Warm Weather Data 

The diurnal effect of solar heating is obvious in 

all warm, weather data. Maximum temperatures are generally 
recorded at 140i0i - 2000 hr, and minima at 020'0 - 0800 hr, with 
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tlie average daily amplitude of variation between 1 and 2®C.. 

Photosynthetic activity should result in maxinium 
DO and maximum pH late in the day (due to the conversion of 
dissolved CO 2 to O2). Plant growth can occur to depths at 
which one percent of the incident sunlight has not been 
absorbed; this is equivalent tO' twice the secchi disk 
reading^ or about 3 to 6 m depth for the St. Clair River. 
Considering this fact, photo^synthetic activity should be 
easily seen for this time period, as the intake was at a 
depth of only 1.5 m. 

This is apparent from the spectral peaks for pH 
in May and DO' in June (Fig. 11) . Less pronounced diurnal 
variations in DO are also evident in April,, with a typical 
amplitude of change being 5 percent in April and up to 25 
percent in June. This magnitude of change in percent satura- 
tion is too large to be explained by temperature changes. 
These observations indicate the effect of rapid weed growth 
on the DO'. 

Variations in the co'nductivity record are less 
regular. In a portion of the April record, maximum conduc- 
tivity and minimum DO tended to complement each other , but 
at other times the time variation of conductivity was irregular. 
These conductivity variations could be due to regular waste 
loadings, causing increased cO'nductivity and oxygen demana, 
but this was not continuously observed. In June, a significant 
sami-diurnal periodicity of conductivity is obvious on, examining 
the spectral plot (Figure 12). However, this periodicity is not 
present in the autocorrelation functio'n (Figure 8) . The 
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reason for this periodicity is not known. 

Calcium and bicarbonate ions are the major co'ntri- 
butors to the conductivity as evidenced by analyses obtained 
by the Great Lakes Surveys (Table 9) . This is generally true 
for most fresh waters (Hutchinson,, 1957, Chapter 9), Data 
for this area of the St. Clair River from the Great Lakes 
Surveys indicate that industrial loadings result in increased 
chloride concentrations near the shore, while midstream data 
are typical of Lake Huron water. The effect of this chloride 
loading on the conductivity is seen in Table 9, The importance 
o,f sulfate ion in determining conductivity is not known in this 
case,, ae analytical data are not available for this ion. 

As the major portion of chloride is due to industrial 
loadings, any periodicity in this parameter will reflect 
periodicities in loading. However, no periodicities are 
observed for chloride at the 80 percent confidence level 
(e,g. Figure 9). If the semi-diurnal conductivity variations 
were due to chloride loadings, a strong semi-diurnal variation 
should have been observed, As this was not the case, the 
source of the conductivity variation is not known. 

4.3.2 Cold Weather Data 

The absence of biological growth during this period 
results in an almoet complete absence of diurnal variations. 
Apparent exceptions are DO in December 1971 and conductivity 
in January 1972. However, the complete spectral plot of DO 
in December indicates a large number o^f peaks below the 80' 
percent confidence level as well as the indicated 24 and 10 
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hour peaks. Interpretations of complex spectra of tliie 

nature must be approached with cautio'n as they may simply 
represent white noise generated by the analyses (Jenkins 
and Watts,, 196 8, p. 157 and 23 0i-9) . The observed periodi- 
cities fO'r pH and conductivity in December can be similarly 
explained. 

In the case O'f January conductivity, a tendency 
towards maximum late in the day and minimum early in the 
day suggests conductivity lo'ading during daylight hours. 
As stated earlier, condoctivity is mainly due to bicarbonate 
but chloride from induetrial wastes is also' imp'Ortant. 
However, no similar trend exists in the chloride data. 
-Sulfate may be important here but analytical data are not 
available. 

Chloride data tend to be random.. For Janijary 
and February 1972, m,ost of the variance is accounted for by 
periods which are close to zerO', where the refined spectral 
estimate values are 1 to- 2 orders of magnitude greater than the 
height O'f the observed peaks. The absolute value of the 
spectral estimate, as well as the size O'f the peak is certainly 
important in determining its significance, and these peaks 
cannot be regarded as significant. 

The 8-hO'Ur periodicity observed for DO in November 
1971 is evident upon examining both the data record and the 
autO'CO'rrelation functio^n. The amplitude of variation is 
frequently 10 percent or more. Some evidence o^f a 4-hO'Ur 
period is alsO' observed, though not at a significance level 
of 8:0 percent. The n^on-normal time distribution O'f DO is alsO' 
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significant, (DO tended to he normally distributed in 

December and January) , This could conceivably suggest 

BOD or COD loading, but other evidence is required. 

Evidence of this nature would be substantiated if data 

were available from several monitors at different lo'catio'ns. 

This would allow a more detailed picture of the waste loadings 

tO' the river . 

In suniiiary^ the periodicities observed during cold 
weather are completely unlike what is seen during the growing 
seasO'n, as is to be expected. They tend to vary from, month 
tO' month, and are accompanied by random short-term changes. 
In many cases, the variance density spectrum is concentrated 
near zero frequencies,^ suggesting that changes are random,,, 
and presumably due to a combination of factors , both natural 
and artificials- 



4.4 Examples of Extreme Data Observed 

One important purpose of automatic water quality 
monitoring equipment is to detect eKtreme results not ordinarily 

detected by routine sam,pling.. This was realized dur,ing the 
spring of 1971 , when evidence suggeeting runoff of salt from 
city streets was obtained. A short period of low dissolved 
oxygen values was also see,n. In the following discussion, 
quoted values refer to edited data as the process of trend 
correction tends to oblite,rate the extremes. 
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The chloride level was highest at the beginning of 
the April record (50 mg/l) and docroasod thereafter with the 
June mean being 23 mg/l. A possible cause o^f the April 
maximum is the discharge of salt used on Sarnia streets to 
the river, follO'Wlng a rainfall. Overflow from the combined 
storm and sanitary sewere may also be important. Meteoro- 
logical observations (Canada Department of Transport, 1971) 
indicate that 0,44 inches of rain fell at Sarnia on April 
12-13 1 two days before the beginning of the record. An 0.29 
inch rainfall on the 17th, however^ was not followed by a 
chloride maxlm^um. An 0.16 inch rainfall on the 2 8th following 
a snowfall on the 24th ^ was followed by a short duration pulse 
of increased chloride on the 30th, suggesting runoff of chloride 
from city streets with a 2-day lag time. Increased chloride 
co'ntent in Toronto Harbour following runo'ff after a snowfall 
has been demonstrated (Brydges , 1969) and increased chloride 
has been evident in Hamilton Harbour fo'llowing spring runoff 
for several years. It is concluded that there is eonsiderable 
evidence for runoff of chloride affecting results at the robot 
monitor, but data are not entirely co^nsistent. 

In May 1971, an unusual decrease in dissolved oxygen 
from 100 to 55 percent saturatio^n was observed o^n May 23-24. 
This was followed by an increase to 7 5-80 percent saturatioin 
on May 25-26, and 100 percent at the end o^f the record (May 27),, 
when calibration indicated no significant error. A comparison 
of edited DO' frequency distributions for April-June 1971, given 
in Figure 13, illustrates the effect of this change on the 
frequency analysis. 
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This is another good c^xample of extreme data 
which would likely be missed during routine surveys. 
However, no correlation with available meteorological 
data existed, and no explanation for this drop in DO is 
available. 



4.5 Markov Chain Analysis 

In order to obtain meaningful results from this 
technique , data selected should have a sufficiently large 
range from minimuiri to maximum that the width of each state 
will be greater than the usual precision of measurement for 
each parameter (Appendix 1, Section Al.4). To this end, 
edited data for chloride, dissolved oxygen and conductivity 
from June and December 1971 were analyzed. The results are 
shown in Figures 14-16. 

The probabilities of occurrence of various states 
for these parameters were determined for each hour up to 40 
hours, given an initial reading within state 1. These prob- 
abilities provide an indication of the hour -by-hour variation 
of these parameters. For example, if the dissolved oxygen 
was 80 percent saturated at a given time (State 1, Figure 15) 
what is the probability of it being 120 percent saturated four 
hours later? The observed probability is close to zero, 
indicating that changes of this magnitude are not observed 
in this length of time. 



24 



In the case of June 1971, figures 14-16 show that, 
following an initial reading in state 1, the state probability 
vector approaches the final state probability vector most rapidly 
for chloride, a good apprO'Ximation being achieved after 12 
hours, (Composite s.ampling would provide the same information). 
However, with conductivity, an extremely slow progression is 
noted, with the probability of being in state 1 still above 
0^. 5 at 4 hours. Even after 24 hours, a good approximation has 
not yet been achieved. This indicates the persistence of 
extreme values of conductivity over periods of even greater 
than one day. If an automatic warning system were established 
to indicate such an extreme, adequate time is available to 
obtain a sample fox further analysis to elucidate the reason 
for such a reading. With dissolved oxygen, a good approxi- 
■mation is achieved after 12 hours, as can be expected from the 
diurnal variation of this parameter . The 4 and 8 hour 
curves indicate some tendency for abnormally low values to 
persist for several hours. 

In, contrast to June 1971, an extremely rapid prO'- 
gression to^ the final state was noted for the December 1971 
conductivity data. (The state vector for 4 hours was a very 
close approximation to the final state vector.) This pre- 
sumably cO'ntrasts the effect of photO' synthesis on the June 
data. On the O'ther hand, the rates of progression for chloride 
and dissolved O'xygen were similar in the two reco'rds, although 
a faster rate of progression might have been expected for 
the December dissolved oxygen, considering the lack of 
photosynthesis in this time period. 



as 
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CODE 81901 


EDITED 








RECORD 


2 














PAGE 2 


DATE 


TINE 


CL 


TEMP(CENT) 


D.O«(f) 


TURB(JCU) 


P7 


PM 


CON0I18 CI CYCLE 


180471, 


14.00 


34.00 


6.89 


112.79 


777777.00 


777777,00 


7,97 


252.96 


180471. 


16.00 


33.00 


6.56 


113.44 


777777.00 


777777.00 


6.30 


254.35 


180471* 


18.00 


31.50 


6,33 


110.88 


777777.00 


777777.00 


8.34 


253.00 


1§*^471. 


20.00 


30.00 


6.17 


110.24 


777777.00 


t77777.00 


8. 21 


252.36 


180471, 


22.10 


32.00 


6.11 


109.44 


777777.00 


777777.00 


8.07 


254.23 


190471, 


O.rt 


33.50 


6.00 


109,60 


777777,00 


777777,00 


8.05 


255.30 


190471 e 


2.00 


2 8. 50 


5.94 


110.08 


7 7 77 77. OO 


777777,00 


8.05 


246.43 


190471. 


4.0O 


30.50 


5.83 


110.08 


777777.00 


777777.00 


8,08 


249,29 


190471. 


6.00 


32.00 


5,9* 


109.60 


777777.00 


777777,00 


8.09 


251.81 


19047lo 


8.00 


28,50 


6.00 


110.88 


777777.00 


777777.00 


8.09 


247, 24 


190471. 


10. on 


28.00 


6.17 


111.68 


777777.00 


777777.00 


8.17 


246.99 


19047lo 


12.00 


29.50 


6.72 


114.59 


777777.00 


777777,00 


8.11 


248.73 


190471. 


14.00 


39.50 


7.06 


114.92 


777777.00 


777777.00 


8,04 


258.07 


190471 o 


16.00 


34.00 


7.00 


115.25 


777777.00 


777777, 0« 


8.11 


255.48 


190471. 


18.00 


31.50 


6,67 


115.57 


777777,00 


777777*00 


3,28 


252.39 


190471. 


20. nf^ 


31.25 


6.28 


113.95 


777777.00 


777777,00 


8.17 


246.61 


190471. 


22.00 


31.00 


6.11 


112.32 


777777.00 


777777,00 


8,07 


247.97 


20n471o 


0.0 


31.00 


6.17 


111.84 


777777.00 


777777,00 


8.08 


248.78 


200471. 


2.00 


31.00 


5.89 


111.68 


777777.00 


777777.00 


8,06 


247,4! 


20f^471. 


4.00 


31.00 


5.83 


112.00 


777777.00 


777777,00 


8,05 


245.71 


20^471, 


6.00 


31.50 


5.89 


111.36 


777777.00 


777777,00 


8,03 


246.52 


200471, 


8.00 


30.00 


6,11 


111.68 


777777.00 


777777.00 


8.06 


246.18 


20O47\^ 


10. OO 


27.00 


6.33 


112.32 


777777.00 


777777.00 


7.96 


240.49 


20047^^ 


1.2.00 


28.00 


6,72 


115.41 


777777.00 


777777.00 


8,06 


244.26 


2004 71a 


14,00 


27.50 


7.C6 


115.57 


777777.00 


777777,00 


8.08 


244.68 


20'^4 71 


16.00 


23.50 


6.83 


116.23 


777777.00 


777777.00 


8.06 


237.84 


2004 71 o 


18,00 


28.50 


6.67 


115,57 


777777.00 


777777,00 


8.23 


249.71 


200471, 


2r*. 00 


33.00 


6, 50 


112.48 


777777.00 


777777. O'O 


8.21 


258.02 


20f^471« 


22.00 


27.50 


6.39 


112.48 


777777.00 


777777,00 


8.12 


245.77 


21 f- 4 71, 


0.0 


30.00 


6,50 


111.68 


777777.00 


777777,00 


8.12 


249.96 


2ir'471. 


2.00 


31.50 


6.39 


111.36 


777777,00 


777777,00 


8,16 


252.92 


2K'47i. 


4.00 


30.00 


6.28 


110.72 


777777.00 


777777,00 


8.19 


249.51 


21^47U 


6.00 


28.50 


6.33 


109,60 


777777.00 


777777,00 


8.20 


246.74 


21C'47U 


8.00 


28.00 


6.44 


109.12 


777777.00 


777777,00 


8,05 


247.48 


21'^47l, 


10,00 


26.00 


6.50 


110.40 


777777.00 


777777.00 


8,01 


245.48 


21.'»471l. 


12.00 


23.50 


6,33 


11<1.56 


777777.00 


777777.00 


8.14 


244.96 



TABLE 1 Sample Robot Monitor Data 





CHLORIOE 




FREQyENCY TABLE 


CLASSFS FREQUENCVIll 
20.00 - 50,50 96, f 3 
50.50 - 31,00 1,75 
91.00 - 111.50 0,0 
111.50 - 142,00 0.0 
14 2,00 - 172.50 O.n 
112,50 - 203,00 ft,0 
201, on - 233.50 0,0 
231,50 - 264,00 0.44 
264.00 - 294,50 0,n 
294.50 - 32 5.00 0,88 


TOTAL 


100, 


f*'O.MENTS 

1 39, 12 '9 
"> 29,36 3 

^ 59, S6? 
4 86,32? 




AVEP^GES FOR THE THiarS 




FISST THIRC = 38.470 




SECOND THIRC = 37,477 




THIRD T'HireO = 30,^54 




THE- HYPOTHFSIS THAT THE S*HFLF IS FRO'l» 
MITH ME4N 13, '3 671 AND VARIAiNCf 
CAN BE REJECTEC HfTH PROB»BILIT¥ 0,0 
' IS '■«5346E 01 FO'P THIS SAMiPLE. 
( 


A INI 'NORMAL D^ISTRI BUTI ON 
945. 3 '9 06 
OF eeiN€ I INCORRECT, THE STATISTIC 1 


ii 

THE HYPOTHESIS THST THE SiMIPLE IS PROP 
'WITH .MiFAN 33,9671 ANO VARIANCE 
CAN BE RfJECTEC WITH PRrieAililTY 0.0 
IS 0« 67691 Ot FOR THIS SA'»1lPLe, 


AIM) EXPONENTIAL Oil ST R I BUT ION 
'945,3906 
CF BEING INC'ORRECT, THE STATISTIC 1 


THE HYPOTHESIS THST THE SAMPLE IS FRO»* 
yiTH ME AM 1.4904 ANlO 'VAiRIANCE 
CAN BE REJECTED WITH PROB.ABIiITY 0,0 
IS 0,3 32 IE 01 FOR THIS SAMPLE. 


A I Ml LCG NORMAL OISTRIBUTION 
0,0184 

OF B.E1NG IMC'ORRECT- THE STATI'STfC Z 



TABLE 2 Raw Data, April 1971 



CHLORIDE 



FRECUEf^Y TABLE 



CLASSES 




20.00 - 


23. 80 


23.80 - 


2 7.60 


27.60 - 


31,40 


31.40 - 


35.20 


35.20 - 


3<?.00 


39.00 - 


42.80 


42.80 - 


46.60 


46.60 - 


50. 4 C 


50.40 - 


54. 20 


54.20 - 


5% 00 



TOTAL 



FBEOUENCYCll 

7.46 
24.12 
37.72 
17.11 

5.26 

3.<»5 

0.88 

1.32 

1.32 

0.88 

100, 



HOMPNTS 




I 


30.43 3 


2 


6.114 


3 


7.253 


4 


9.921 



o 



AVERAGES FOR THE THIRDS 
FIRST THIRD = 34.648 

SECOND THIRD ^ 28.727 

THIRD THIRD ^ 28.197 



THE HYPOTHESIS THAT THE SAMPLE IS FROM AINI NORMAL DISTRIBUTION 
WITH MEAN 30.5241 AND ¥ARIANCF 34.2983 

CAN BE REJECTEC yiTH PROBABILITY O„000 OF BEING INCORRECT. THE STATISTIC Z 

1% 0„ 2604E 01 FO'R THIS SAMPLE. 

THE HYPOTHESIS THST THE SAMPLE IS FROM A(NI EXPONENTIAL DISTRIBUTION 
WITH MEAN 30.5241 AND VARIANCE 34.2983 

CAN BE REJECTEO WITH PROBABILITY 0.OO2 OF BEING INCORRECT. THE STATISTIC Z 
IS 0.1854E 01 FOR THIS SAMPLE. 



THE HYPOTHESIS THST THE SAMPLE IS FROM A(N) LCG NORMAL DISTRIBUTION 
WITH MEAN 1.4774 AND VARIANCE 0.0059 

CAN BE REJECTED WITH PROBABILITY 0.001 OF BEING INCORRECT. THE STATISTIC Z 
IS 0.1974E 01 FOR THIS SAMPLE. 



TABLE 3 Edited Data, April 1971 





cm LOR IDE 






FiiQUENCY TAiiE 


' 


CLASSES 


FREQUiNCYfll 


iuii - 


23.69 


1.12 


5»3.69 - 


25.67 


5. 26 


25,61' - 


2 7,65 


11,84 


27.65 - 


29.63 


21.49 


29.63 - 


31.61 


22.81 


31.61, - 


33. 5 i 


2 2. 81 


33,58 - 


35.56 


7, B9 


35.56 - 


37,54 


3, 51 


3T,54 - 


39.5? 


1,32 


3^.52 - 


41. 5ft 


1,75 ; 


TOTAL 




lOfi, 


MOMENTS 






\ 30.5 57 






1 3,35? 






1 2.352 




, 


* 4,574 






AVERAGeS FO'R THE THIRDS 






FIRST THfPC ^ ■il1.4!6 






SECOND THI^C ^ 30, I2fl 






THiRDi THIPO ^ 30.177 






THE HYPaiTHESrS TH4T THE SMPLE IS FROhM' 


i 1 IN 1 NO R *l A L 1 S TR I ByT 1 ON 


WITH mm 30,. §401 ,ANC VARIAIC'P 


11.3331 




CAN Bc PiJECTEC WITH PROJAB^liL IT¥ 0i,673 


OF SEING IMCORRECT. 


THE STATISTIC Z 


IS O'.iatiE m for this sam^ple. 






THE HfPOTiNESIS THST THE ,SAMPLE IS FR01M 


A! Ml EXPOiNElNlTIAi OrSTR liyTION 


WITH iM'FftM 1fi.S401 AMD' VARIANCE 


11.3331 




CAN iE BEJECTEC WITH PRQiBABJiLITY €.Q.OO 


CF SEIFiG INCORRECT, 


THE STATISTIC I 


1% ^.IflTE O'l FOR THIS SAHiPie. 






THE HYPOTHESIS THST THE SAMPLE IS FRC'N 


kim LOG NORHAi DISTRISUTIONI 1 


: WITH mAH 1.4122 ANiO WARIANCE 


0.0023 


"=■ 


CAN iBE REJECTED WITH PROBABILITY 0.T52 


OF iilMG INCORRECT, 


THE STATISTIC I 


IS o*675 3E 00 FOR THIS SAMPLE. 







TABLE 4 Trend Corrected Data, April 1971 



TABiiE 5 
POLYNOMIAL FT'TTING DURING TREND CORRECTION 



Data Period 


Number of 
Intervals 


Level 


of . 


Bignlficance 
(percent) 


from F-Test 








Cl° 


Temp 


m 


DO 


PH 


Cond . 


April 1971 


i 




1.0 


NS 




1.0 


NS 


NS 


May 1971 


f 




- 


5.0 




* 


NS 


NS 


June 1971 


13 




ili. 


0.1 




5.0 


mm 


5.0 



early Nov. 1971 

December 1971 

January 197 2 

February 197.2 



11 



0.1 



NS 



2.5 



14 


NS 


0.5 


0.1 


NS 


_ 


10 


NS 


2.5 


0.5 


NS 


NS 


i 


* 


■*■ 


NS 


NS 


NS 


13 


1.0 


NS 


_ 


NS 


2.5 



Note : * Although the pO'lynomlal was not significant at 5%,, 

trend correction was accepted because the portion of 

the variance attributed tO' deviation about regres- 
sio^n was lees than 2% O'f the total variance. 

NS indicates that the po'lynomial fit was not significant 

at the 5% level. 

The chloride probe was out of calibration in May and 

early November 1971. The pH reference probe was mal- 
functioning in June 1971 and calibration in November 
was questiO'nable. The dissolved oxygen probe was out 
of calibration in February 1972. 
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TABLE 6 

MEANS AND STANDARD DEVIATIONS OF CORUNNA 
ROBOT MONITOR DATA 



Data Period 

April 1971 

May 1971 

June 1971 

early 
November 1971 

December 1971 

January 1972 

February 1972 



CI (mg/l) 
Mean S . D 



30.5 



22.8 



3.7 



6.1 



32.0 6.0 

29.1 8.9 
35.0 8.7 



Temp. (°C) 
Mean S . D . 



6.3 0.4 

10.7 0.7 
13,0 0.6 

12.8 1.0 
6.2 1.1 
1.6 0.8 
1 . 5 . 7 



Dissolved Oxygen 
(% saturation) 
Mean S,D. 


pH 
Mean 


S.D. 


Conductivi 

Ciimho/cm at 1 

Mean S . 


105.7 


2.2 


8,2 


0.1 


253 


8 


83.1 


10.1 


8.1 


0.3 


317 


9 


103,1 


5.8 


*t^' 


•m, 


310 


7 


82.3 


5..g 


mm 


fa 


224 


9 


84.5 


4.9 


8.1 


0.3 


235 


13 


74.0 


4.8 


7.7 


. 2 


228 


22 


mm 


lw-< 


7.8 


0.3 


241 


14 



Notej_ Trend removed data are used where significant, otherwise edited data. 

Readings were obtained at hourly intervals, except in April 1971, when 
the interval was 2 hours. 



TABLE 7 
DISTRIBUTION PROBABILITIES BY KOLMOGOROV-SMIRNOV TEST 



Data 
Period 



April 1971 



May 1971 



June 1971 



early Nov 
1971 

Dec 1971 



Jan 1972 



Feb 197 2 



CI 



Temp . 



DO 



PH 



Conductivity 



Dist. Prob. Dist. Prob. Dist. Prob. Dist. Prob. Dist. Prob 



log ,7 52 
norm. 



log .0'12 
norni. 



log ,055 normal .22 9 
no^rm . 

normal .0 00 log .172 

norm. 

normal .632 normal .000 



log .002 normal .094 
norm. 



normal .02 5 normal .0 



no'rmal .790 log .527 normal .8 36 

norm , 

noxmal .04 5 normal .000 



log .306 exp, 
norm. 



noxmal .0 



normal .24 6 normal .0 



normal . 



.016 



log .184 
norm . 

normal .000 



normal .453 



log .994 
norm. 



normal .002 log .8 30 

norm. 



normal .405 



Note ; Trend r'emoved data has been used where significant, otherwise edited, 
^"~^~ The distribution with the highest probability is indicated. 



TABLE 8 
80% CONFIDENCE SPECTRAL MAXIMA 



No. of 
Data No, of 



Data Period Points Lags Cl Temp. DO pH Cond. 



April 1971 228 30 None *24 hr (24 hr) (24, (24, 

12 hr 12 hr) 12 hr) 

May 1971 162 30 - 24 hr None 24 hr None 

June 1971 625 100 None *24 hr *24 hr - 12 hr 

early Nov. 525 30 - None 8 hr - None 
1971 

Dec 1971 689 100 None None 24 hr 10 hr 16 hr 

10 hr 6.3 hr 

Jan 1972 144 30 5 hr None None None 24 hr. 

Feb 197 2 900 30 10 hr None - None None 



Notes : * indicates that the peak is significant at the 95% 

confidence level. 

( ) indicates interesting periodicities which were 
visible but below the 80% level. 

Data were prewhitened by binomial weighted running 
mean over 5 hourly observations. 
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TABLE 9 

SOME ANALYTICAL DATA FOR ST. CLAIR RIVER 
DETERMINED BY GREAT LAKES SURVEYS 



This data is for Station SR-31.0, 42-53-54N, 82-28-18W, 
1 Km north of the robot monitor intake. 









Conductivi 


^ty 
















(25°C) 




CI 


Alkalinity 


Date 




Location 


pmhos 




mg/1 


mg 


CaCO; 


3/1 


May 


5/71 


M 
S 


213 
274 




7 
26 




•tmr 




May 


31/71 


M 
S 


209 
256 




5 
19 




81 
92 




June 


2 8/71 


M 
S 


212 
252 




6 
17 




89 

88 




Oct. 


26/71 


M 


20 9 




7 




84 








S 


276 




24 




88 





Note : M = midstream (640 m from U.S. shore) 

S = Canadian shore (1^13 m from U.S. shore) 

Average equivalent concentrations: CI (M) 0.2 meq/1 

(S) 0.6 meq/1 
Alk 1 , 7 meq/1 
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PC P 




LOCATION OF 

ROBOT IWONITOR INTAKE 



iCKWilTH Sf" 



M^L 



HILL 



ST 



CORUNNA 



INTAKE AT 

LAT 42^53' 25" N 

LON a 2® 2 7' 20" W 



Fig, 1 Location of robot monitor 
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i 
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Fig. 2 Raw chloride^. April 1971, frequency distribution 
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Fig. 3 Edited chloride, April 1971, frequency distribution 
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Fig.. 5 Autocorrelatio^n coef f icients foT June 1971 chloride 
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Fig. 6 Autocorrelation coefficients for June 1971 temperature 



1.0 |- 



0.8 ~ 



ae 



-0.2 - 



ul 






m 






u. 






y. 




X 


O 0.4 


— 




« 






m 






«. 




X 


■1^. 






< 






_i 






^ 0.2 


— 


X 


QC 






O 






O 






© 






■|- 




X 


=> 






< 






0.0 




X « 
X 
^ X 






^xx^ 



x^'^x 

X X 

X X 



vXX 

X* X 

X X 



X X 

X ^x 
xxx 



XXXv 

X "x 
X X 

1 X 

«. -1 

X X 

X' X 

X V ' 

„ » X 

X X XX 

XX X X 

XvX y X 



xXx 

X x 



0.4 >— 





J \ 1 L 



10 20 



30 



40 50 60 

LAG TIME -HOURS 



70 80 90 too 



Fig. 7 Autocorrelation coefficients for June 1971 dissolved oxygen 
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Fig. 8 Autocorrelation coefficients for June 1971 conductivity 
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Fig. 9 Variance density spectra for June 1971 chloride 
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Fig. 10 Variance density spectra for June 1971 t^ein,peratu 



re 



46 



,1 r- 



.05 



jUIi; 

w 

m 



>• 

m 

m 

a 

H 
O 






.01 



.005 






.001 



.0005 



.0001 



XX 



I 

\ 
I 
I 

X 

1 
i 

I 

4 



I 



V,' 

AY 



x^ 
I t 

I I 

X 

I 

I 



95%-' '-90% 



V 



X 






Si}i 



I 

i 
\ 



x^x 



^x 

X \ 

I 



I 

\ XH 

I I . 

I I t 
X / X 



0.05 



0.10 O'.IS 

FREQUENCY tir^l 



0.20 



0.25 



Fig. 11 Variance density spectra for June 1971 dissolved O'Xygen 
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Fig- 12 Variance density spectra for June 1971 conductivity 
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Fig. 13 Dissolved oxygen frequency distributions for April-June 1971 
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Fig. 14 State probability distributions for chloride,. June 1971 
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APPENDIX 1 

DESCRIPTION OF THE ROBOT MONITOR 
AMD PROBLEMS ENCOUNTERED DURING ITS OPERATION 



hl.l Description of Installation 

The "robot monitor" unit is located in the 
Corunna waste treatment plant. The intake and pumping 
arrangement, as of August 1971, is shown in Figure Al.l. 

The intake line is 30 meters of 3 cm diameter 
pvc hose^ connected to 4 m of 10 cm diameter perforated 
pvc pipe directed downstreaiti^ and anchored to a steel 
tripod at a depth of 3.5 m. The piece of large tubing 
served to avoid uptake of weeds, sand or other debris on 
the river bottom. Concrete blocks were also attached along 
the hose at approximately 2m intervals, to keep the hose on 
the river bottom. This helps to prevent dragging of the line 
due to river currents resulting from passing ships or ice 
floes. A pump of 1 1/2 h.p. is used to provide a constant 
flow of 8-12 liters/min required for each cell. This pump, 
with an electric heater to maintain an adequate indoor temp- 
erature, is housed in a wooden hut sitting on a culvert. 
As this unit is entirely in the water, it is subject to 
action of currents and waves, as well as ice in the winter, 
and a more permanent concrete housing has been considered. 

A 2.5 cm steel pipe abo^ut 4 m, long connects the 
pump to a 2.5 cm plastic hose about 80 m long, buried at a 
depth of 0,5 - 1 m in porous soil, leading to the robot monitor. 
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Although freezing prO'blems might be expected in the winter, 
difficulties were encountered only once (February 1971) in 
the fO'Ur years the installation has been in use. 



A1.2 Problems Associated with Installation 

(1) Weeds 

On July 29, 1971, the original intake line, which 
was about 10' m lo^ng at a depth of 1,5 m, was observed to be 

in a concentration of weed growth. The growth was worse on 
August 16,, when the line was extended past the weeds into an 
area of clear water. However, on August 23,, the flow had 
dropped from 20' to 6 1/min due to weeds plugging the intake 
line. At this time, the 4 m length O'f perforated 10 cm, 
diameter pvc tubing was added, and flow conditions became 
normal . 

(2) Ice 

Action of ice in February and March 1972 resulted 
in the intake line being broken and carried away. Sand was 
taken in,, and at one point the DO and pH cells were plugged. 
Before the extent o^f damage was known,, attempts were made to 
keep the system in operation by backf lushing once or twice 
a week. On May 29,, the line was replaced, and proper flow 
conditions were again achieved. 
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Al . 3 Problems Associated with the Instrument 

(1) Digitizer 

The digitizer was newly installed on April 14 , 
1971, Problems which have been encountered are as follows 

(a) At times the alignment of the punch head 

was incorrect. Frequently , the punch holes 
were at the edge O'f the tape. Because O'f 
this ,, a replacement punch head was used in 
the period May-July 1971, while the o-riginal 
was repaired. Afterward, occasional adjust- 
ment was required, 

(b) The timer controaiing the take-up motor is 

not reliable. This no,rmally allows the take-up 
motox tO' run for a set interval following 
punching O'f a set of data, in order to wind 
up the punched tape. Failures o^f this unit 
have resulted in continuous operation of the 
take-up mo'tor (Sept. 2 6/71) and failure of 
the take-up motor to operate (June 24, 1972), 
In bO'th cases, the tlm,er was replaced. 
Presently, this problem, is being discussed 
with Schneider Instrument Co. 
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(2) He liability of Probes 

NO' problems were encountered with the temperature 
and conductivity prO'bes, with the changes at calibratiO'n 
being within the expected precisions. 

The DO probe presented considerable problems. 
As a result of silver chloride build-up o^n the silver anodes, 
which flaked off from time to time during operation, as well 
as slime deposition on the membrane, frequent reconditioning 
and recalibration was required. Comparison o^f digitizer 
reading with results obtained by Winkler titration indicated 
errors o^f + 10 percent or more in the period April 1971- 
January 1972, Large negative errors were consistently 
found in February -March 1972, but this occurred when the 
intake was malfunctioning. In order to improve the reliability, 
the calibratiO'n of the probe should be checked weekly against 
the Winkler results, as was done after March 3, 1972. Weekly 
reports by S, Bechard before this date only give panel meter 
readings; it is suggested that digitizer readings, as well 
as chemical test results be reported weekly for all parameters. 

The chloride probe was reliable over most of the 
time period, with only small adjustments being needed on re-^' 
calibration. A nO'table exception was in the May 4-May 27, 
1971 periO'd, in which the ratiO' digitizer. analyzer reading 
was 3 to 3.5, rather than the expected 2. It was stated that 
calibration was perfO'rmed but nO' reason for the observed dif- 
ferences was given. 
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Reliability of the pH probe was comparable tO' 
the chloride probe. At times, crystalLizatiO'n of KCl in 
the reference electrode resulted in unreliable results such 
as in June 1971. A serious error also was noted, on November 
22, 1971, when the pH 8 buffer read 9.10'. At other times, 
recalibration adjustments amounted to + 0.1 - 0.2 pH units. 



Al.4 Accuracy of Equipment 

Estimated ace^uxacies, as claimed by Schneider 
Instrument Company, are as follows: temperature + 1°F; 

cO'nductivity , DO and pH, + 1%; chloride,, + 5%. The ranges 
of measurement were: cO'nductivity , 0-1200 and 0-2400 
micromhos; DO, 0-12 and 0-24 mg/l; chloride, 0-600 mg/1. 



Al . 5 Maintenance of Equipment 

Periodic visits to Corunna have been performed 
by J. Gabay, and the results of these visits recorded in the 

log book,,. During these visits, all sensors and other equip- 
ment were checked, and required service performed as suggested, 
in the i,nstruction manual. ,Meter readings of the various 

parameters were dompa,red, with digitizer readings, and the 
results reccrded . Where necessai-y, recalib,rat,ions were per- 
f orm,ed . From Septeiiibe,r 17, 1971, weekly checks were iTia,de 
by S. Bee hard. 
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APPENDIX 2 
COMMEHT ON QUALITY OF VARIOUS DATA RECORDS 

TO' facilitate diecussiO'n of the best data O'btained, 
from the Corunna robot monitor^ connments on all data obtained 
have been included in Table A2 . 1 , The best reco'rds were 
obtained fox' April-June 1971 and December 1971, and these 
are discussed in detail elsewhere, along with any data from 
several "fair" reco^rds , as applicable. In July 1971, cali- 
bration difficultieB were experienced. In August and 
September 1971^ the paper tape was punched off center, 
although the September tape was translated following re- 
punching. The November 1971 record was only fair due to 
calibration difficulties. 

During 1972,^ very little, If any, good data were 
obtained. Problems with ice, which have been mentioned in 
Appendix 1, resulted in approximately three months' data 
being obtained from the stream bottom, where considerable 
mud was drawn in by the pump, clogging sensors several times. 
The data for February 1972, which were obtained prior to the" 
start of ice difficulties , are believed to be good, except 
for dissolved oxygen, which was out of calibration. During 
the period May tO' July 1972, electronic difficulties of gradu- 
ally increasing magnitude were experienced. Failure of a 
stepping relay in the digitizer during this time period 
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resulted in records of alternating high and low readings, 

which could not be translated meaningfully. At the end 
of July 1972, the instrument was taken out of service for 
repairs. 
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Location 
Code 



81901 

81903 

81904 

81905 
81906 
89107 

8190 8 

81909 

81910 

81911 

81912 
81913 
81914 

ti * 

81915 
81916 

81917 

81918 



TABLE A2.1 



INDEX O'F DATA TAPES OBTAINED AT CORU'NNA 



Period of 
Operation 



Text 
Reference 



Comments 



Apr 15 ~ May 4/71 



April 1971 



May 20 - May 2 7/71 May 1971 

May 28 - June 23/71 June 1971 

June 23 - July 28/71 July 1971 

July 2 9 - Aug 16/71 August 1971 

Aug 17 '- Sept 14/71 Sept 1971 



Sept 30 = Oct 14/71 
Oct 14 - Oct 2 5/71 
Oct 2 5 - Nov 16/71 
Nov 16 - Nov 26/71 

Nov 26 - Dec 25/71 
Jan 14 - Jan 20/7 2 
Jan 20 = Feb 27/7 2 



early Oct, 
1971 

late Oct. 
1971 

early Nov. 
1971 

late NO'V. 
1971 



Dec 1971 

Jan 197 2 

Feb 197 2 



Mar 2 - Apr 7/7 2 Mar 197 2 

Apr 7 - Apr 17/72 Apr 1972 

Apr 17 ~ May 30/72 May 1972 

May 30 - June 25/72 June 1972 

June 29 - July 27/72 July 1972 



good record, 2 hr. 
frequency 

good record r Cl out 
of calib. 

good record , except 

PH 

Cl~r DO ^ pH no good 

tape punched O'ff center 

fair recoxd ,, no good 
periodicities 

poor record, high 

deviations 

poor record f high 
deviations 

pH and Cl~ out of 
calibration 

Temp and DO question- 
able. pH out of 
calibration 

gO'Od reco'rd 

fair record 

fair record, DO no 
good 

intake at bottom 

intake at bO'ttO'in 

intake at bottom 

poo^r record , instru- 
ment needed overhaul 

as for June 1972 



Note: The text reference is an approximate date reference used in 
— ~~^ referring to the various tapes in the text and other tables. 

* This record was split into two portions due to its extreme 
length. 
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APPENDIX 3 
STATISTICAL ANALYSIS OF FREQUENCY DISTRIBUTIONS 

A comparison of several statistical methods of 
analysis of robot monitor data for normal , exponential and 
log normal distribution was performed. The following methods 

we re s e 1 e c ted . 

1« Kolmogorov-Smirnov (continuous distribwtion) 
2.. Von Mises (equation for no^rmal distribution) 

3. Von Mises (discrete distribution) 

4. x^ test (discrete distribution) 

KOLMOGOROV-SMIRNOV TEST 

The Kolmogorov-Smirnov test, which is designed for 

comparing twO' distribution functions, is the test which has 

previously been used (IBM,. 196 8) . In this test, the data 
(expressed in terms of number of standard deviations from the 
mean) are sorted into ascending order. The maximum deviation 

D from the theoretical distribution is then obtained,, .and the 
test statistic Z is calculated as : 

Z = D /nT [M] 

where n = number O'f readings. 

To test the effect of sample size on the probability 
a,nd to compare this test with the Von Mises test fo,r normal 
distribution, a set of chlo'ride data was sampled and tested for 
normal a.nd exponential distribu,tion. The results are give,n i,n 
Table A3.1. 
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The value of Z for exponential distribution is 
approxiniately inversely proportional to the square root of 
tlie sampling frequency. For normal dietribution, the rate 
of decrease is less. At any rate, for very small numbers of 
readings, the probability of normal distribution also increases 
rapdily. This illustrates the hazard of attempting distri- 
bution calculations on small data samples. 

VON MISES TEST (Equation for Normal Distribution) 

Von Mises (1964) devised a test which is applicable 
to both discrete and continuous distributions. The observed 
distribution function S (x) is compared to the theoretical 
distribution function P(x), and e is defined as i 

e (x) = n [S(x) - P(x)) [A2] 

where n is the number of observations. As a measure of the 
deviation, a parameter u^ was defined as: 

if? = A(x) e^(x)dx [A3] 

where X is a weight function, supposedly evaluated to make the 

expected value of oj- equal to 1. In the case of the normal 

distribution, it can be shown (Von Mises, 1964) that: 

n 
= 2 /? I ( <l)(x ) + x^ [<I>{x^) -(2i-l)/2n] ] =n [A4] 

i=l 



2 
CD 



where $ (x) = ^ / ^"^^ dx, the equation for the 

/2Y '"'°° 2 no'rmal distribution 

and <p (k) = ' (x) 
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The standard deviation O'f w^ is given by: 

o = /0,63 - 0.12/n [A51 

Finally, the probability of finding the normal distribution is 
calculated from the normal probability curve for the deviation 

O'f u)- {in term,s of a) from its expected value. 

The Von Mises test for normal distribiition was 
checked for several samples. It was found, however, that 

E(a)^) = 0, rather than 1, as stated by Vo^n Mises. This was 
verified by taking a data set of 100 readings of normal 
distribution, taken from the normal probability table. For 

this data lo^ = 0.0078. Some typical probabilities are lower 
than those obtained by the Kolmogorov-Smirnov test, and,, also 
increase with decreasing sample size. As well, this method 

is not suited tO' calculating probabilities of exponential 
distribution. Therefore, its use was not further considered:. 



VO'N M,I;SES TEST FOR A DISCRETE DISTRIBUTION 

The data are divided into a given number of classes, 
each class representing an equal fractio^n of the difference 
between the minimum and maxim.um reading. The number of obser^ 



vations n. within each class i is determined. u^ is then calcu- 
1 

lated from the following equation. 

k (nQ. - N, ) 
i=l 



2 _ .., 11" . /^F [A6] 
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where N. = En. = the number of observations in 

j=l ^ all classes below and including 
class i 

Q. - the normal cumulative distribution function 

fO'r the mean of class i 

B ^ the total number of readings 

a ^ the standard deviation of all readings 

k s^ the total number of classes. 

The standard deviation of w^ is calculated from equation [A5] 
as in the continuous case, and a probability of normal distri- 
bution can be estimated from the normal curve. If desired, 
the probability of exponential distribution may be estimated 

by fitting Q,, from an expo^nential curve. 

X 

Edited chloride data from December 1971 were subjected 
to this test. The results are given in Table A3 . 2 . The value 
of 0)^ is directly proportional to the same size n, as ie 
evident from the equation. AlsO', w^ increases as the number 
of classes into which the sample is divided is increased. These 
drawbacks produce severe limitations on the validity of using 
the Von Mises test for discrete distributions on the present data 
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CHI -SQUARE TEST (Discrete Distribution) 

A frequency distribution is obtained as discussed 
for the Von Mises test. A, nDxirial distributiO'n is fitted, 

and expected frequencies obtained. The value of x^ is then 
obtained from the equation (Snedecor and Cochran, 1967): 

k 
X' = 2 (f. - F. )2 IA71 

.=1 -^— 

where k - number of classes 

f. = observed frequency in class i 
.F,. = calculated frequency in class i. 

The probability is calculated from the x^' distribution with 

(n-S) degrees of freedom. Fitting the exponential distri- 
bution presents the drawback that terms in x^ become infinite 
for data more than one standard deviation below the mean. 
This difficulty was circumvented by replacing F.. in the denomi- 
nator of the above equation by f . when F. was close to or equal 
to zero . 

Resulte O'f this test on edited chloride data from 
December 1971 are presented in Table A3. 2... As in the Von, Mises 
test, the value of x^ is directly proportional to the number 
of observations J the effect o^n the probability is even wo^rse 
than that observed with the Vo^n Mises test, and allows the 
probability tO' assume values in almO'St the complete range of 
to 1. The value of x^ also increases with increasing number 
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of classes, but this change is compensated for by the increased 
number of degrees of freedom* 

As a result of these drawbacks, it was decided to 
continue using the Kolmogorov-Smirnov test for estimating 
distribution probabilities . 
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SUMMARY 

The performance of several statistical tests for 
evaluating probabilities of normal, exponential and lo-g- 
no^rmal distribution, was evaluated. All methods suffer from 

the drawback that the obs^erved prob^abilities are dependent 

upon the sample size, but this dependence is worst for the 
discrete distribution tests. The following additio'nal draw- 
backs are no'ted : 

1. In the Von M,ises tests (both continuous and 
discrete), the expected value of w^ is 0, not 
1 as claimed by Von Ml,ses. 

2. The Von Mises continuou,s test is unsuitable 
for exponential distribution. 

3. Probabilities for the discrete Von Mises test 
decrease as the number O'f classes in the 
frequency table is increased. 

As the KolmogO'rov-Sm,irnov test is least affected by 
these drawbacks, it will continue tO' be used in estimating the 
distribution probabilities . 
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TAB.LE_ A3 . 1 

KOLMOGORQV-SMIRNQV AND VON MISES TEST 
RESULTS ON SAMPLED DATA 

D.ata used are edited CI from January 1972., (144 readings) 



Sampling 
Frequency 




Exponential 
Z PrO'b 


No 
Z 


rmal 
Prob 


Von 


Mises 
Prob 


[E(a)'2)^ 


-0] 


% 




2 


.083 


.0003 


.850 


.466 


1.041 


.18 9 






t 




1 


.482 


.025 


.593 


.874 


.4 56 


.565 






1 




1 


.167 


.131 


,561 


.911 


.4 02 


.612 






5 






.929 


.355 


.4 33 


.992 


,245 


.757 






■i 






.817 


.518 


.420 


.995 


.229 


.7 73 






10 






.66 8 


.763 


.369 


.999 


.161 


.838 






12 






.577 


.8 93 


.336 


• 9999 


.132 


.8 67 







if 



TABLE A3. 2 



DISCRETE NORMAL DISTRIBUTION TESTS 



Data used are edited CI from Deceniber 1971 (689 readings) 



*4 



No. of 
Classes 

20 
20 
10 
10 



Method of 

Expressing 

Frequencies 

No . observations 

percent 

No . observations 

percent 



Von Mises 
co^ Prob [E(tja^)^0] 

1.114 0.159 

0.162 0.838 

0.294 0.710 

0.0427 0.957 



Chi-Square 



xl 

46.24 
6.71 

19.91 
2.89 



Prob 

0.0002 
0.987 
0.006 
0.895 



Note : The Kolmogorov-Smirnov probability for nonnal distribution 
of the above data is 0.025. 
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APPENDIX 4 
OPERATIONS AND CALIBRATIONS LOG 



1, Operatione Log 



Date 



Service 



April 14/71 

May 19/71 

May 2 7/71 
June 23/71 

July 16/71 
July 28/71 

Aug. 16/71 

Aug, 2 5/71 
Sept. 15/71 

Sept. 26/71 

Sept. 30/71 
Nov. 22/71 

Jan. 13/72 
Mar. 6/72 



Digitizer inetalledi DO probe 

replaced . 

Digitizer skipping punches and 

off center; alignment corrected. 

Digitizer punch head replaced. 

DO probe replaced. 

pH and CI reference electrodes 
filled with KCl solution. 

DO probe replaced. 

Original punch head replaced. 
Intake in weed growth. 

Intake extended, anchored with 

cement blocks . 

pH and Cl~ reference electrodes 

cleaned, filled. 

Intake extended, 10 cm pvc pipe 
of length 3.5 m installed. 

pH and Cl~ reference electrode 

cleaned, filled. 
DO probe replaced. 

Digitizer malfunctions: improper 
spacing of punched tape, takeup 
motor running continuously. 
Timer replaced. 

Punch unit replaced, spacing adjusted 
Conductivity sensor repaired, 

DO probe replaced . 

pH and CI- reference electrodes 

cleaned and filled. 

Adjusted punch unit, eerviced 

digitizer. 

Fixed open ground on pH probe 

DO probe replaced. 

Intake line removed by ice. 

Frequent, backf lushing of intake 

performed. 
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APPENDIX (continued) 



April 17/72 

May 29/72 

June 29/72 
Aug. 4/72 



Attempt to repair intake 
unsuccessful due tO' ice 

movement in river . 

Stepping relay cleaned^ 

intake line repaired. 

pH and CI"" reference electrodes 
cleaned and filled. 

Timer of takeup motor replaced. 

Stepping relay cleaned. 

MO'nitor taken out o^f service 
for repairs. 



f-«' Sample Calibration Log. (May 27/71) 

Note I RO'Utine calibration was performed O'n each of the 
above dates. 



Parameter 



Analyzer Digitizer 
Reading Reading 



Measurements 
Taken 



Temperature 
Conductivity 

'BO 

CI" 

pH 



Sl'^'F 512 
355 pmho 354 



12.75 mg/1 599 (=11.98 

mg/1) 

30' mg/1 10'6 (-53 

mg/1) 



9.6 



838 



52.7«'F 

350i pmho; calibrated 

on standards 
140, 680, 
970 pmho 

12.8 mgyi 



50 rog/1 



Calibrated on buffers 
pH=7 and 4 



?1 



